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Abstract.  The  most  recent  solution  of  Earth  orientation  parameters 
from  the  observations  made  by  optical  astrometry  in  1899.7-1992.0  at  33 
observatories  is  used  to  estimate  the  corrections  of  the  present  IAU  model 
of  precession-nutation.  Since  the  resolution  of  the  solution  is  about  five 
days,  only  the  terms  with  periods  of  14  days  and  longer  are  considered. 
The  results  are  compared  with  VLBI-based  corrections  of  the  standard 
precession-nutation  model,  and  the  differences  are  discussed. 


1.  Introduction 

Earth  Orientation  Parameters  (EOP  -  polar  motion  x,y,  Universal  Time  UT1- 
TAI  and  celestial  pole  offsets  Ae,A ip)  have  been  recently  derived  from  optical 
astrometry  observations  made  since  the  beginning  of  the  century  at  33  observato¬ 
ries.  The  Hipparcos  Catalog  (Kovalevsky  et  ah,  1997)  was  used  to  recalculate  all 
the  observations;  they  were  thus  brought  into  the  recently  adopted  International 
Celestial  Reference  System  ICRS  (Feissel  &  Mignard,  1998).  The  complete  de¬ 
scription  of  the  algorithms  applied  in  our  preceding  solution  are  described  by 
Vondrak  et  al.  (1998a).  Small  changes  introduced  later  to  produce  the  most  re¬ 
cent  solution  are  given  in  Vondrak  et  al.  (2000).  The  changes  between  these  two 
solutions  consist  in  including  data  from  two  more  observatories  (Mount  Stromlo 
and  Jozefoslaw),  adding  more  data  from  Kharkov  and  Blagoveschtchensk,  dis¬ 
carding  the  data  from  Ukiah  after  1960,  and  more  Hipparcos  proper  motions 
corrected  than  before.  These  solutions  are  here  denoted  as  OA97  and  OA99. 


2.  The  solution 

The  details  of  the  solution  OA99  being  described  in  detail  elsewhere,  only  a  short 
description  is  given  here.  The  solution,  covering  the  interval  1899.7-1992.0,  is 
based  on  approximately  4.5  million  individual  values  of  three  types  of  observ¬ 
ables  (instantaneous  latitude  ip,  difference  of  Universal  Time  from  the  uniform 
atomic  time  scale  (UT0-TAI)  and  the  difference  between  observed  and  calcu¬ 
lated  zenith  distance  of  the  star  8h),  according  to  the  type  of  instrument  used. 
The  observations  were  made  with  50  different  instruments  at  33  observatories. 
In  case  two  or  more  similar  instruments  were  used  at  the  same  observatory, 
their  results  were  merged  into  a  single  series,  with  the  steps  due  to  different  lo¬ 
cations  of  the  instruments  removed,  and  treated  further  as  a  single  instrument. 
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Figure  1.  Celestial  pole  offsets  from  optical  astrometry  OA99. 


In  addition,  we  derived  also  a  set  of  auxiliary  ‘station’  parameters,  giving  the 
systematic  deviations  of  individual  instruments  and  the  rheological  parameter 
A  =  1  +  k  —  l  (responsible  for  nonrigid  tidal  variations  of  local  verticals).  The 
list  of  the  instruments  used  in  the  solution  is  given  in  Table  1. 


3.  Celestial  pole  offsets  and  their  analysis 

Because  we  used  the  standard  IAII  models  of  precession  IAU1976  (Lieske  et 
al .,  1977)  and  nutation  IAU1980  (Seidelmann,  1982)  to  calculate  the  apparent 
positions  of  the  observed  stars,  the  celestial  pole  offsets  derived  from  our  solution 
express  the  displacements  of  the  Earth’s  pole  from  the  celestial  epliemeris  pole 
whose  position  is  defined  by  these  models.  The  results  are  plotted  in  Figure  1, 
the  values  themselves  in  the  upper  part  (crosses),  their  formal  uncertainties 
below  (lines),  both  in  the  same  scale.  It  is  necessary  to  say  that  A £{,  A 
referred  to  the  same  epoch  are  strongly  correlated  (their  error  ellipses  highly 
elongated),  due  to  only  a  partial  coverage  of  observations  during  a  day  (for  a 
detailed  explanation  of  this  effect  see  Vondrak  et  al.  (1992).  The  error  ellipses 
are  always  pointed  towards  the  Sun  (the  observations  being  centred  around 
local  midnight)  and  they  thus  rotate  in  the  celestial  reference  frame  with  an 
annual  period  that  causes  an  annual  change  in  rms  errors  of  As  and  At/’  (phase 
shifted  by  90°).  This  effect  became  much  weaker  after  1956,  when  latitude 
observations  (the  only  type  of  observations  used  before  that  date)  were  combined 
with  Universal  Time  observations.  The  figure  shows  also  that  the  noise  level 
substantially  diminishes  after  1956,  thanks  to  a  larger  number  of  participating 
observatories  (and,  consequently,  also  observations). 

The  series  of  celestial  pole  offsets  was  subject  to  analysis,  in  which  we  used 
a  classical  least-squares  method  to  estimate  corrections  to  precession  plus  the 
nine  most  significant  nutation  terms  (OA99/1).  Since  the  resolution  of  the  series 
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Table  1.  List  of  participating  observatories  and  instruments  in  OA99. 


Observatory 

Instr. 

Interval 

long. 

lat. 

Beijing 

PAST 

1979-87 

116°20' 

40°06' 

Belgrade 

ZT 

1949-85 

20  31 

44  48 

Blagoveschtchensk 

ZT 

1959-91 

127  30 

50  19 

Bratislava 

CZ 

1987-91 

17  07 

48  09 

Carloforte 

ZT 

1899-43,46-78 

8  19 

39  08 

Cincinnati 

ZT 

1899-16 

84  25 

39  08 

Gaithersburg 

ZT 

1899-14,32-78 

77  12 

39  08 

Grasse 

PAST 

1983-91 

6  56 

43  45 

Irkutsk 

ZT,PTI 

1958-90,79-91 

104  20 

52  17 

Jozefoslaw 

ZT 

1961-91 

21  00 

52  06 

Kharkov 

PTI 

1973-91 

36  14 

50  00 

Kitab 

ZT 

1930-78 

66  53 

39  08 

Mizusawa 

ZT,FZT 

PZT#1,2 

1900-78,67-84 

1959-75,74-91 

141  08 

55 

39  08 

55 

Mount  Stromlo 

PZT 

1957-85 

149  00 

-35  19 

Nikolaiev 

PTI 

1974-91 

31  59 

46  58 

Ondrejov 

PZT 

1973-91 

14  47 

49  55 

Paris 

AST 

1956-82 

2  20 

48  50 

Pecny 

CZ 

1970-91 

14  47 

49  55 

Poltava 

ZT#1,2,3 

1949-90,50-68,68-80 

34  30 

49  36 

Praha 

CZ 

1980-84,85-91 

14  25 

50  05 

Pulkovo 

ZT 

PTI#1,2,3 

1904-41,48-91 

1959-71,71-85,71-91 

30  20 

55 

59  46 

55 

Punta  Indio 

PZT 

1971-84 

57  17 

35  21 

Richmond 

PZT#1,2 

1949-87,81-89 

80  23 

25  36 

Santiago 

AST 

1965-90 

70  33 

33  24 

Shaanxi 

PAST#  1,2 

1974-84,85-91 

109  33 

34  57 

Shanghai 

AST, PAST 

1962-84,75-84 

121  26 

31  11 

Simeiz 

AST 

1977-90 

34  00 

44  24 

Tschardjui 

ZT 

1899-19 

63  29 

39  08 

Tuorla- Turku 

VZT 

1963-89 

22  30 

60  25 

Ukiah 

ZT 

1899-60 

123  13 

39  08 

Washington 

PZT#1,2,3 

1915-55,54-84,81-91 

77  04 

38  55 

Wuchang 

AST, PTI 

1964-86,81-86 

114  21 

30  32 

Yunnan 

PAST 

1980-91 

102  48 

25  02 

ZT  -  visual  zenith-telescope,  VZT  -  visual  zenith  tube,  FZT  -  floating  zenith-telescope, 
PZT  -  photographic  zenith  tube,  PTI  -  photoelectric  transit  instrument, 

AST  -  Danjon  astrolabe,  PAST  -  photoelectric  astrolabe,  CZ  -  circumzenithal 
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is  live  days,  only  the  long-periodic  terms  (between  13.6  days  and  18.6  years)  are 
considered.  In  order  to  estimate  how  robust  the  solution  is,  we  did  the  same 
analysis  with  the  most  precise  part  of  the  solution,  covering  only  1956.0-1992.0 
(OA99/2).  Both  sets  are  further  compared  with  the  values  obtained  by  analyzing 
a  recent  VLBI  solution  of  GSFC  (covering  only  a  twenty-year  interval,  1979.59- 
1999.92),  using  exactly  the  same  algorithm.  The  results  of  all  three  analyses  are 
depicted  in  Table  2.  The  correlations  between  individual  values  (except  for  the 
constant  term  that  is  always  correlated  with  the  trend  because  of  the  choice  of 
the  epoch  J2000.0)  in  the  table  are  rather  small  for  optical  astrometry  (They  do 
not  exceed  0.13  in  case  of  OA99/1  and  0.29  in  case  of  OA99/2.),  thanks  to  the 
length  of  the  series.  It  is  not  the  case  for  VLBI  where  the  correlation  is  as  high 
as  0.95  between  the  trend  and  cosine  term  of  the  18-yr  period. 


Table  2.  Corrections  to  precession  and  long-periodic  nutation  terms, 
derived  from  optical  astrometry  observations  (OA99/1  from  1899.7- 
1992.0,  OA99/2  from  1956.0-1992.0)  and  VLBI.  In  the  first  three  rows, 
sine  columns  give  the  trends  in  milliarcseconds  per  year,  cosine  columns 
the  mean  celestial  offsets  at  epoch  J2000.0. 


Argument  / 
Period  [days] 

Longitude 
sin  (‘in’) 

[mas] 

cos  (‘out’) 

Obliquity 
cos  (‘in’) 

[mas] 
sin  (‘out’) 

0 

oo 

OA99/1 

OA99/2 

VLBI 

—  2.156±0.046 
— 1.516±0.143 

—  2.990±0.013 

—  31.02il.71 

—  16.15i3.52 
-43.00i0.10 

-9.17i0.64  - 

-7.12il.38 
— 4.90i0.05  - 

■0.093i0.018 

0.007i0.056 

■0.220i0.007 

Q 

6798 

OA99/1 

OA99/2 

VLBI 

-4.31±1.36 
-3.84il.72 
—  6.89i0.03 

-2.15il.39 
—  2.26il.81 
3.39i0.07 

1.83i0.50 

1.71i0.69 

2.89i0.04 

2.76i0.50 

2.36i0.66 

1.44i0.01 

2Q 

3399 

OA99/1 

OA99/2 

VLBI 

3.12il.32 

2.28il.65 

1.08i0.03 

2.41il.33 
2.97il.67 
—  0.15i0.02 

-1.10i0.49 
-1.34i0.66 
—  0.29i0.01 

1.22i0.49 

1.51i0.65 

0.10i0.01 

2F-2D+2Q, 

182.6 

OA99/1 

OA99/2 

VLBI 

-14.50il.36 

-13.16il.68 

1.57i0.02 

-4.03il.39 
-4.12il.72 
—  1.36i0.02 

4.68i0.50 
4.45i0.66 
—  0.58i0.01 

—  2.28i0.50 
-0.82i0.65 

—  0.43i0.01 

v 

365.3 

OA99/1 

OA99/2 

VLBI 

-14.15il.23 

-4.82il.54 

4.99i0.01 

6.58il.56 

8.50il.90 

1.24i0.01 

10.69i0.48 

9.04i0.63 

2.10i0.01 

5.99i0.52 
6.19i0.67 
—  0.27i0.01 

V+2F-2D+2Q, 

121.7 

OA99/1 

OA99/2 

VLBI 

-3.82il.36 

-2.75il.69 

O.OliO.Ol 

-0.46il.37 

—  1.93il.70 

—  0.02i0.01 

-1.92i0.49 

-0.54i0.65 

0.04i0.01 

— 0.08i0.50 
0.82i0.65 
—  0.05i0.01 

2F-2D+Q, 

177.8 

OA99/1 

OA99/2 

VLBI 

3.31il.36 
3.76il.70 
—  0.13i0.02 

2.97il.37 

4.45il.70 

0. OliO. 02 

— 3.69i0.50 
-3.80i0.65 
O.lliO.Ol 

2.80i0.50 

2.75i0.65 

O.OliO.Ol 

2T+2Q 

13.66 

OA99/1 

OA99/2 

VLBI 

-0.18il.33 
-0.90il.66 
—  0.35i0.02 

1.37il.32 

1.36il.66 

0.29i0.02 

0.58i0.49 

0.68i0.65 

0.14i0.01 

—  072i0.49 
-0.99i0.65 
0.17i0.01 

l 

27.55 

OA99/1 

OA99/2 

VLBI 

-0.09il.32 
-0.92il.65 
—  0.13i0.01 

-0.47il.32 
-0.25il.65 
—  0.06i0.01 

-0.73i0.49 

-0.69i0.65 

0.06i0.01 

1.04i0.49 

1.20i0.65 

0.06i0.01 

2  F+tl 

13.63 

OA99/1 

OA99/2 

VLBI 

2.23il.33 
3.15il.67 
—  0.19i0.02 

-l.91il.32 
-l.19il.65 
—  0.02i0.02 

0.46i0.49 

0.91i0.65 

0.07i0.01 

—  0.14i0.50 
-0.23i0.65 
0.05i0.01 
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4.  Discussion  and  conclusions 

First  of  all,  one  should  notice  rather  big  differences  among  the  trends  (sine 
columns,  first  three  rows)  in  longitude:  2.156(OA99/l),  1.516(OA99/2)  and 
2.990(VLBI)  and  in  obliquity:  -0.093(OA99/1),  0.007(OA99/2)  and  0.220(VLBI) 
These  differences  are  definitely  much  larger  than  their  formal  errors  but,  on  the 
other  hand,  the  differences  between  OA99/1  and  VLBI  are  smaller  than  those  of 
our  preceding  solution  OA97  (Vondrak  et  al. ,  1998b)  that  yielded  1.543  mas/yr 
in  longitude  and  0.131  mas/yr  in  obliquity.  The  correction  to  precession  is  now 
(solution  OA99/1)  much  closer  to  the  VLBI  value  than  before.  This  change  is 
very  probably  due  to  our  less  tolerant  approach  to  Hipparcos  proper  motions 
(20%  of  them  being  corrected  in  OA99  in  contrast  to  only  10%  in  OA97).  Also 
interesting  is  the  obvious  change  of  trend  in  longitude  around  1956,  noticed  al¬ 
ready  in  the  OA97  solution  by  Bizouard  et  al.  (1997).  This  effect  can  probably 
be  also  caused  by  imperfections  of  the  Hipparcos  proper  motions.  (Different 
stars  were  observed  in  the  second  half  of  the  century  than  in  the  first  one.) 

Secondly,  there  is  a  discrepancy  in  the  cosine  (out-of-phase)  term  of  the 
principal  nutation  in  longitude  between  both  optical  solutions  and  VLBI  (they 
differ  in  sign).  This  is  the  very  term  that  is  highly  correlated  with  the  trend 
in  the  VLBI  solution.  This  also  holds  for  the  cosine  term  in  longitude  with 
a  period  of  3398  days  (see  also  Yaya  et  al.,  (2000)  for  a  similar  discussion). 
The  largest  differences  between  optical  astrometry  and  VLBI  are  however  found 
in  the  semiannual  and  annual  terms.  In  addition,  the  values  found  from  the 
present  solution  OA99  differ  substantially  from  OA97.  It  is  clear  that  this  is  the 
weakest  point  of  optical  astrometry.  We  suspect  that  this  effect  is  mostly  due  to 
proper  motions  of  some  of  the  stars.  Different  stars  at  different  right  ascensions 
are  observed  during  the  year  (Only  the  ‘night’  stars  are  observed.),  and  if  their 
positions  differ  from  reality  at  a  certain  part  of  the  sky  they  can  produce  an 
annual  wave  in  celestial  pole  offsets.  All  the  other  terms  agree  quite  well,  on  the 
level  of  their  formal  standard  errors. 
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